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Abstract. - We construct a new spin-1 model on a chain. Its ground state is determined 
exactly which is three-fold degenerate by breaking translational invariance. Thus we have 
trimerization. Excited states cannot be obtained exactly, but we determine a few low-lying ones 
by using trial states, among them solitons. 



> 

in 

m 

o 

a^. 

^^ ' Spin chains and spin ladders have been intensively studied in recent years |l|. Besides 

the Bethe ansatz solvable models there is an increasing number of models where at least the 

ground state is known exactly and where eventually the gap to the first excited state can be 

calculated Q . Typically in the ground state of these models the correlations are of short range 

^-^ • which implies a finite gap. One natural way to construct such models in one dimension is to 

^ I use the matrix product method [Q which is based on the concept of " optimum ground states" 

O ■ as explained in ref . |j] . 

In some cases translational symmetry is spontaneously broken leading to twofold degenerate, 
dimerized ground states. To the best knowledge of the authors no model has been found so 
K^ , far where the ground state would be threefold degenerate, thus leading to trimerization. For 

$^ ' the first time we construct such a model. 

There exists one model, the Uimin-Lai-Sutherland (ULS) model M, which is a special case 
of the spin-1 bilinear- biquadratic chain, where the excitation spectrum has a tripled periodicity 
in the Brillouin zone, namely soft modes appear at q = and q = ±27r/3, however, the ground 
state itself does not break translational symmetry. 

For later convenience the Hamiltonian of the spin-1 bilinear-biquadratic model is written 
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in the form 

N 



n^Y. [(^'^ • s»+i)' + " (S'' • S'+i)] • (1) 




The ULS point corresponds to a = 1. For a > 1 a gap opens up in the excitation spectrum 
above the unique ground state, if the chain is periodically closed. At a = 3 we get the AKLT 
model 0, where the existence of the gap has been proven. For a < 1 the gapless spectrum with 
three soft modes presumably survives M, 0]. Since a gapless spectrum and, correspondingly, 
power law like decaying correlations should be an exception, and generically a spin chain should 
have a finite gap, one might expect to find a gapped trimerized state by perturbing away from 
the Hamiltonian given above. 

To construct such a model that has an exact threefold degenerate ground state with a finite 
gap to the excitations we observe that the quantity 

A,,, = (S, • S,)' + a (S, • Sj) +a-l (2) 

has the following properties. As eigenstates arrange in multiplets, we denote the 2-spin singlet 
state at site i and j by s[z, j], and similarly by t[i,i] and q[i,j] the triplet and quintuplet states. 
Then we find 

(3-a)s[i,j] , 

Oi[i,j], (3) 

2aq[i,j]. 

This implies that the triplet has the lowest energy in the range < a < 3. 

For three spins at sites i,j, fc one has one septuplet, two quintuplets, three triplets, and 
precisely one totally antisymmetric singlet. This trimer singlet is given by 

5[z,j,fc]--^[(+,0,-) + (0,-,+) + (-,+,0)-(-,0,+)-(0,+,-)-(+,-,0)] (4) 

in terms of the three states (+), (— ) and (0) of a spin-1 operator. Then using the fact that 
in order to form a singlet with the third spin, any two spins have to form a triplet, one gets 
easily from eq. o) 

A J S[i, j, k] = A^^k S[i,j, k] = Aj^k S[i, j, k] = 0. (5) 

Using this property one sees immediately that the Hamiltonian 

N N 

'H — 2_^ /i'i, 1+1,4+2,1+3 ~ 2_^ Ai^i+2Ai+i^i+z (6) 

i=l i=l 

has precisely three zero energy ground states. For a finite periodically closed chain of iV = 3p 
sites, where p is an integer, these three states are: 

p 
*i = J]^ S-pi - 2, 3i, 3i + 2] , 

i=l 
P 

*2 = Y[S[3i~l,3i + l,3i + 3], (7) 

i=l 
P 

$3 = Y[S[3i,3i + 2,31 + 4]. 
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They are ground states, as Ti. and h of eq. (|6|) are positive-semidefinite in < a < 3, and they 
are optimum ground states, as in the sense of ref. [^ they are simuhaneously ground states of 
aU local interactions h. For a finite system these states are not strictly orthogonal, but their 
overlap goes to zero exponentially with N ^ oo namely as 7^ (0 < 7 < 1). 

To visualize the states we show one of them in fig. m by connecting the sites belonging to 
the trimer singlets by valence bonds. Because the trimer singlet is antisymmetric under left 
to right reflection, i.e., [i,j,k] — > [fc,j, z], the bonds have to be directed. It is more convenient 
to draw the chain as a zig-zag ladder which is shown in fig. 0. When N = 6p, this is a genuine 
ladder, while when N = 3(2p +1), it can be thought of as part of a Moebius ribbon. In what 
follows we will use everywhere this zig-zag ladder representation. 



Fig. 1. - Valence bond representation of the ground state. 



Fig. 2. - The ground state configuration represented on a zig-zag ladder. 



That these states are in fact the ground states of the Hamiltonian in the range < a < 3, 
follows immediately from eq. 0. At a = the quintuplet of the spins at site i and j has 
the same energy as the singlet, so at this point the degeneracy of the ground state becomes 
exponentially large. For a < the quintuplet has the lowest energy, leading to the fully 
polarized ferromagnetic state. The transition of the trimer phase to the ferromagnetic one is 
of first order. 

At a = 3 the 2-spin singlet and triplet local pairs become degenerate, producing again an 
exponentially large degeneracy, while for a > 3 we get to the massive, Haldane-like regime 
with unique ground state. 

Next we discuss excited states of the model. They cannot be calculated exactly, but clearly 
there has to be a gap to the excitations, because the ground state is formed of short range 
valence bonds. The simplest excited state could be a state where one trimer singlet is broken 
up and is replaced by a triplet (Tri), quintuplet (Quint) or septuplet (Sept) of the three spins, 
e.g., 



$b(0 = B[il- 2,31,31 + 2] YlS[3i-2,3i,3i + 2], 



(8) 



where B stands for Tri, Quint or Sept. Using these as trial states and calculating their energy 
as the average E = (Ti.) , we see that two terms of the Hamiltonian give a non-zero contribution: 



^3/-2,3;^3/-i,3;+i + ^3;-i,3/+i^3;,3;+2 



(9) 



As can be easily seen the Hamiltonian has no matrix element between states where different 
trimers are broken, and there is no overlap between these states. Therefore these excited states 
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have no dispersion. For the three possible triplet states we get the energy 

r (2 + a)(i + a), 
£^Tri = < (3-ai -Q!)(i +a), 
[ (3 — 02 — a)(| + a) , 

where 01^2 are the solutions of the equation 



(10) 



ai,2 = 1- §a± a/(1 
For the two quintuplet states the energies are 



i-y 



-E, 



Quint 



3a(i + a) . 



(U) 



(12) 



while for the septuplct state Egcpt = 4a(l/3 + a). One checks that the triplet gap with ai in 
eq. (|l^) and the first quintuplet gap in eq. ( |l^ are the lowest ones. 

The lowest energy propagating modes are probably solitonic excitations, which are moving 
domain walls between two different ground states. In a ring of length N = 3p, where the 
lattice spacing is taken to be unity, a soliton cannot appear alone. In the simplest case one 
needs three solitons to satisfy the periodic boundary condition. However, provided the chain 
is long enough, solitons can be studied separately assuming that the two ends are in different 
ground states. 

A domain wall that perturbes the system as little as possible can be obtained on a chain of 
3p + 2 sites between state ^1 on the left hand side and V&s on the right hand side by having 
the spins at sites 3/ — 2 and 3^ + 1 forming a singlet, triplet or quintuplet, while all other spins 
are in their trimer singlet ground state. Similar domain walls appear between state ^2 on the 
left and ^1 on the right or ^3 on the left and 'I'2 on the right. 

Such a state can be written as 



;-i 



$c(0 = n '5'[3« - 2, 3i, 3i + 2] c[3l - 2, 3/ + 1] J| S[3j, 3j + 2, 3j + 4] , 



(13) 



i=l 



j=l 



where c stands for s, t or q. It is shown in fig. ^. 



31+1 




Fig. 3. - Domain wall on a chain with 3p + 2 sites represented on a zig-zag ladder. 



Only one term in the Hamiltonian will give a non-zero contribution to the energy, namely 

^3/-2,3;^3/-l,3; + l • (14) 

Although the Hamiltonian has no matrix element between states with different site index Z, a 
propagating mode is obtained since there is a finite overlap between these states. Now looking 
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for the moving soliton in the form 



'i>c(g) = j]e^^''$,(?) 



1=1 



we find for the energy of the singlet and quintuplet excitations 



EsAq) 



4 \3 



[5 — 3cos3q], 



and 



for the triplet excitations. 



^t('?) = i(i+")'[5 + 3cos3q] 



(15) 

(16) 
(17) 



31+1 




31-2 



31+4 



Fig. 4. - Larger domain wall on a chain with 3p + 2 sites represented on a zig-zag ladder. 



A somewhat more complicated domain wall is shown in fig. m between state ^i on the left 
and state ^3 on the right. Because of its symmetric shape, the Hamiltonian has finite matrix 
elements between domain walls on different sites, giving direct propagation. In addition to that 
the overlap has also to be considered, giving finally a rather complicated dispersion relation. 
When the two spins form a singlet, the dispersion relation of the soliton is: 



Es{q) 



(l + 3a) +(l-2a) cos3(7 



5 — 3 cos 3q 
14/3 -2 cos 3g' 



while for the triplet 



^t(<z) = 2(i+a)'[l + icos39] 



5 + 3 cos 3q 
16/3 + 4 cos 3g' 



and for the quintuplet 

E^{q) 



(l + 3a)Vi(l + a)^cos3<7 



5 — 3 cos 3g 
14/3 -2 cos 3<?' 



(18) 



(19) 



(20) 



It is interesting to speculate what happens if the model presented in this paper is perturbed 
away from the form given in eq. m). As mentioned earlier both at a = and a = 3 the 
degeneracy of the ground state gets exponentially large. Beyond these points the system 
becomes ferromagnetic, or a first order phase transition to a massive Haldane phase occurs. If 
extra parameters are introduced by making the couplings of the two-spin and four-spin terms 
independent, the trimerized phase will survive because of its finite gap. However, at some 
point the two legs of the zig-zag ladder become decoupled, and other phases may appear. 

Finally we mention that the model can easily be generalized without changing essential 
properties. Namely for N — 6p one can have different parameters a and (3 in the operators A 
acting on the upper and the lower legs, respectively. We like to thank Dr. A. Schadschneider 
for this hint. 
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Furthermore we mention that we have found other model Hamiltonians with trimer ground 
states. However, in all these cases they are degenerate with other ground states and the 
degeneracy is exponentially large. Clearly they are not so interesting. 
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